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ABSTRACT 

We present deep 1.4 GHz Very Large Array (VLA) radio continuum observations 
of two ~half square degree fields in the Coma cluster of galaxies. The two fields, 
"Coma 1" and "Coma 3," correspond to the cluster core and southwest infall region 
and were selected on account of abundant pre-existing multiwavelength data. In their 
most sensitive regions the radio data reach 22 /iJy rms per 4'.'4 beam, sufficient to detect 
(at 5a) Coma member galaxies with L1.4 ghz = 1-3 x 10^'' W Hz~^. The full catalog 
of radio detections is presented herein and consists of 1030 sources detected at > 5a, 
628 of which are within the combined Coma 1 and Coma 3 area. We also provide 
optical identifications of the radio sources using data from the Sloan Digital Sky Survey 
(SDSS). The depth of the radio observations allows us to detect AGN in cluster elliptical 
galaxies with Mr < —20.5 (AB magnitudes), including radio detections for all cluster 
ellipticals with Mr < —21.8. At fainter optical magnitudes (—20.5 < Mr < —19) the 
radio sources are associated with star-forming galaxies with star formation rates as low 
as 0.1 Mq yr-^ 

Subject headings: galaxies: clusters: general — galaxies: clusters: individual (Abell 
1656) — galaxies: radio continuum 



Introduction 



The Coma cluster of galaxies, Abell 1656 (jAbell. Corwin. Qlowinlll989l . hereafter re ferred to 
as Coma), has been among the most well-studied galaxy clusters (e.g., see the review by iBiviano 
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19981 ). At z = 0.0231 (jStruble &: Roodlll999l ) it is nearby, yet significantly richer than Virgo and 
more accessible than Perseus (A426), which lies near the Galactic Plane. Studies of Coma have 
shaped our understanding of galaxy evolution in the cluster environment. In assessing the frequency 
of blue galaxies in intermediate redshift clusters, Coma wa s one of several nearby cl usters used as a 
benchmark for establishing the "Butcher-Oemler Effect" (iButcher &: Oemlerlll984l ). Although the 
fraction of blue galaxies decreases with decreasing redshift and therefore the activity level of Coma 
is pr e sumab ly lo w, studies have st il l revealed g alaxy evolution in progress. ICaldwell et al.l (jl993l . 



1996 



19991 ) and ICaldwell &: Rosd (|l997l . Il998l ) have studied images and spectroscopy of Coma 
galaxies to reveal populations of starbursts and post-starburst galaxies, and related these findings 
to substructures within Coma and hence connected the star-formation history of galaxies to cluster 
environment. In particular. Coma member starburst and post-starburst galaxies a,re preferentially 
locat e d in the regions aroun d NGC 4839 which appears to be a large infalling group (jColless &: Dunn 



1996 



Neumann et al.l 120031 ) althou gh others have ar gued for this system being seen shortly after 



passage through the core of Coma (jBurns et al.lll994l l. Similar studies have contrasted Coma with 
higher redshift clusters; at higher redshift mor e luminous galaxies are observed to constitute the 
starburst and post-starburst populations (e.g., iDressler et al.lll999l). while source s in Coma with 
similar spectroscopic signatures are fainter dwarf galaxies (jPoggianti et al.ll2004l ). This suggests 
that "downsizing" - the shift in star-formation activity to lower mass systems - is occurring, and 
that some of the Butcher-Oemler effect is the result of applying a fixed absolute magnitude cutoff 
to both low and high redshift systems. 

Rec ently, Coma has been studied to greater depth using the William Hersch el Telescope 
(WIIT ). lKomiyama et al.l (120021 ) presented optical imaging in five 32f5x50f8 fields, with lMobasher et al 
(|200ll ) using the imaging to select targets for spectroscopy for the "core" and "outskirt" fields: 
"Coma 1" centered at 12'^59""45!2 -F27°57'53'.'l, and "Coma 3" at 12'^57'^28!5 -F27°08'08'.'5 and 
having NGC 4839 near its northern edge. The southwest corner of the Coma 1 field overlaps slightly 
with the northeast corner of the Coma 3 field, meaning that combined they define a contiguous 
total area of about 0.92 square degrees (refer to Figure [1]). These p hotometric and spectr oscopic 
data were used to construct the cluster o ptical luminosity funct ion (jMobasher et al. 2003), study 
stellar population ages and metallicities ( Poggianti et al. 2001al fbl). investigate the e ffect of env i- 



romnent on galaxy propert ies and compare the findings to higher redshift clusters ([Carter et al 



2002 




Poeeianti et al. 


(Edwards et al. 


2OO2I) 



Poggianti et al.l l2004l ) , and examine cluster dynamics for both giant galaxies and dwarfs 



Surprisingly, the best wide area 1.4 GHz radio continuum su rvey of Coma to date is the Fain t 
Images of the Radio Sky at Twenty centimeters survey (FIRST: iBecker. White. Sz Helfandlll995l ). 
With a limiting sensitivity of 1 mJy and 5" resolution, the FIRST data are capable of detecting 
Coma member sources with Li,4 ghz ^ 1-2 x 10^^ W Hz~^. This amounts to radio luminosities 



about half that of the Milky Way and greater (|Condonlll992l ). or in terms of a star formation rate 
(SFR) about 0.7 M(n yr~^ (as suming a Salpeter IMF from 0.1 to 100 Mq, using the relationship of 



Yun. Reddy. &: Condon 



2001 



While this is a useful limit for studies of active galaxy populations 
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in clusters (e.g.. iMiller &: Qwenll2002l . identify an average of '^10 radio-detected galaxies per cluster 
down to this limit and within each cluster's Abell radius), the proximity of Coma makes it relatively 
easy to study hitherto unexplored fainter radio galaxy populations. Deeper 1.4 GHz observations 
in the vicini ty of NGC 4 8 69, t he prototypical head-tail radio galaxy at the core of Coma, are 
presented in iFeretti et al.l (|l995l ). These data reach an rms sensitivity of about 30 /iJy at about 
the same resolution as FIRST, and have the added benefit of full polarization information. They 
are of very limited areal coverage (a few square arcminutes), however, on account of bandwidth 
smearing. Similarly, the excellent H i observations of spiral galax ies in Coma have narrow total 



bandwidths and hence relatively poor radio continuum sensitivity (jBravo-Alfaro et al.l 120001 . 12001 



the rms for 1.4 GHz continuum of these data is about 0.2 mJy). These data show evidence for cluster 
environmental effects in the form of ram pressure stripping, including asymmetric H i morphologies 
and greater gas deficiencies for galaxies near the core of the cluster. 

There is also a wealth of new multiwavelength data availa ble for Coma. X-ray satellites have 
imaged both the core and large-scale regioiis of the cluster (e.g., 



White et al. 



1993 



Briel et al. 



2001 



Vikhlinin et al. 2001 : Neumann et al. 2003 : Finoguenov et al. 20041), and much deeper observa tions 



aimed at the Coma 3 region have been performed with Chandra (jHornschemeier et al.ll2006l ) and 
XMM (Hornschemeier et al., in preparation). Fewer prior data exist in th e ultraviolet (UV ) , 
although Coma was observed by the FOG A balloon-borne imaging telescope (jPonas et al.lll995l ). 
However, new UV data from GALEX including a deep observation centered on Coma 3 are now 
available and iiidicate a faint-end slope to the UV LF that is steeper than comparable field UV LFs 
(jHammer et al.ll2008l ). Although more traditionally associated with star formation, the UV data 
from GALEX detect numerous elliptical galaxies, even down to less luminous, dwarf-like objects. 
The reason fo r this is emission from ex treme horizontal branch stars, known as the "UV upturn" in 
their spectra (jCreggio &: Renzinilll99d ). In additi on to the aforemen tioned WHT optical programs, 
a new Hubble Space Telescope Treasury program (jCarter et al.ll2008l ) obtained two-color imaging of 
portions of Corn a using the Advanced Camera for Surveys. The Sloan Digital Sky Survey (SDSS; 
York et al.ll2000l ) include s optical imaging and spectroscopy fo r Coma in its fifth and sixth data 
releases (DR5 and DR6; lAdelman-McCarthv et al.l boOTI . 



20081 1 



In the infrared (IR), Spitzer is 
producing i nteresting resu l ts on the cluster. Near-IR observations with IRAC (3.6 - 8.0 /xm) are 
reported in iJenkins et al.l (|2007l ) for the Coma 1 and Coma 3 regions, and also suggest a steep 
faint-end slope for the near-IR LF. These LF results are indicati ve of an excess o f dwarf galaxies 
relative to optically-selected samples. At longer IR wavelengths, iBai et al.l (|2006l ) constructed 24 
/im and 70 fxm LFs using Spitzer MIPS and found their overall shape to be consistent with field 
LFs. They did, however, find evidence for environmental effects in the form of different slopes for 
the core and infall regions. 

To complement these deep and wide area multi- waveband datasets, we used the National Radio 
Astronomy Observatory (NRAofl Very Large Array (VLA) to obtain deep 1.4 GHz images of the 
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core (Coma 1) and outskirts (Coma 3) of the Coma cluster. Our observations reach 22 /iJy, or 
LiA GHz ~ 1-3 X 10^'' W Hz~^ for galaxies at the dista nce of Coma. For Co ma members powered 
by star formation, this equates to SFR > 0.08 Mq yr~^. iHopkins et al.l ( 20021 ) find the median SFR 



of blue compact dwarf galaxies as determined by 1.4 GHz emission to be 0.3 Mq yr~^, implying 
that our VLA survey should be sufficiently sensitive to provide radio detections for some of the 
brighter dwarfs, and useful upper limits for non-detections. The radio source catalog derived from 
our VLA observations has been correlated with the SDSS optical data to explore the radio galaxy 
population of Coma, and by extension numerous radio sources lying behind the cluster. 

In this paper, we present the observations, data reduction, and resulting source catalogs. A re- 
lated paper investigates the characteri stics of the radio-d etected Coma member galaxies and derives 
the cluster radio luminosity function (jMiller et al.ll2008l ). For all distance-dependent calculations, 
we assume a redshift of z = 0.0231 corresp onding to a distance modulus of 35 using the stan- 



dard WMAP cosmology (jSpergel et al 



2007). This is consistent with distances to Coma derived 



using a host of othe r techniques (e.g., ICreggl EqQtI ; Hjorth Sz Tanvir 1997 : Kavelaars et al. 
Tullv fc Pier^l200nl ;i. 
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2. Radio Data 
2.1. Observations 

The VLA observations were performed over five days in June 2006 as program code AM868. 
On each of the five days the scheduled time was centered on the transit of Coma. The VLA was 
in its B, or second largest, configuration. In order to obtain nearly-uniform sensitivity across the 
~ 30' X ~ 50' areas corresponding to Coma 1 and Coma 3, mosaics of individual pointings were 
made. These amounted to 11 separate pointings for each field, arranged in a hexagonal grid with 
12'.5 between adjacent pointings (see Tableland Figure [I]). This grid spacing insures that all 
regions of the Coma 1 and Coma 3 fields with the exception of small areas in six of the corners 
lie within the half-power response of more than one pointing. Observations of 3C286 for primary 
flux density calibration and 1310+323 for bandpass and phase calibration meant that each of the 
22 total pointings received just under one hour of total integration. To minimize the time spent 
slewing while still obtaining good (u,v) coverage, we devoted individual days to Coma 1 and Coma 
3 (2 days each) and reobserved all 22 pointings on the fifth day. On days when a single field was 
targetted we observed each of its 11 pointings twice for about 11.5 minutes each time, with the 
order shuffled on the other date that field was observed. Thus, in total each pointing received five 
individual ~ 11.5 minute observations at five different hour angles, resulting in excellent total (u,v) 
coverage as demonstrated in Figure [2] for the centermost pointing within the Coma 1 field. 
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The correlator was operated in the standard multichannel continuum mode used for deep, 
wide-field 1.4 GHz observations. This corresponds to seven 3.125 MHz channels in each of two 
bands (called "intermediate frequencies" or IFs) centered at 1.365 and 1.435 GHz, respectively, 
with each IF receiving both left and right circular polarizations (hence 28 total channels). The 
smaller channels alleviate bandwidth smearing while still providing a fairly large total bandwidth 
and thus increasing the sensitivity. At the time of the observations the VLA had just begun its 
transition to the Expanded VLA ("EVLA"), with individual VLA antennas being retrofitted with 
EVLA receivers and electronics. Depending on the date of the observations either three or four 
EVLA anntennas were used in the observations on a shared-risk basis, along with 22 standard 
VLA antennas. The EVLA antennas were operated in the same manner as the standard VLA 
antennas (e.g., same observed frequency channels), so the benefits of including them were simply 
the increased number of array elements. 



2.2. Data Reduction 

The data were reduced using the Astronomical Image Processing Software (AIPS) package, 
version 31DEC06. The flux density scale was tied to 3C286, which was bootstrapped to the phase 
calibrator 1310-1-323. This latter calibrator was observed approximately every 45 minutes and also 
was used for calibrating the bandpasses. Some flagging of data was necessary, usually for the EVLA 
antennas. For example, one EVLA antenna (the one most recently brought into service) provided 
useable data for only one IF and polarization. The others were reliable with the exception of short 
periods of time usually at the conclusion of each day's observations. 

Once the data were calibrated, the initial steps toward producing final images were made by 
pointing ID and observation date. Thus, there were 66 separate (u,v) datasets with which to work, 
corresponding to 11 pointings each for the first four observing dates and all 22 pointings on the 
last (refer to Table [1]). The data for each pointing were imaged using the task IMAGR with a 
strong taper to produce an image 2° on a side, onto which was overlaid a "flys-eye" pattern of 37 
slightly overlapping 512" x 512" square facets. The apparent positions of bright sources outside the 
coverage of the flys-eye were no ted, and their true positions were taken from the FIRST catalog 



(jBecker. White. Sz Helfandlll995l ). There were usually around ten such bright sources located out- 
side the central flys-eye pattern. In subsequent imaging, all 37 facets of the flys-eye and all the 
bright sources outside of it were imaged separately using the 3D option of IMAGR. Additionally, a 
few of the brightest sources within the flys-eye coverage were also given their own facets for imaging. 

Each of the 66 (u,v) datasets were then subjected to iterative rounds of imaging and self 
calibration. After initially imaging the ~ 50 facets for a specific pointing, the imaged facets were 
inspected and sources were boxed for subsequent imaging rounds. The imaged facets were used 
as inputs for self calibration and additional editting of the (u,v) data. The end product was 66 
well-behaved (u,v) datasets, and the three separate days of data per pointing were then combined 
to form the semi-final 22 (u,v) datasets (i.e., one per mosaic pointing). Each of these were again 
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subjected to imaging and self calibration, with the final images per pointing having sensitivities 
ranging from 39 to 45 /iJy per 4'.'4 beam. The observational strategy ensured nearly identical (u,v) 
coverage for each pointing, resulting in circular 4'.'4 restoring beams for all pointings. 

We also performed quick tests of the relative performance of the retrofitted EVLA antennas. 
These were achieved by imaging the data for a given pointing both with and without the EVLA 
antennas included. The rms noise of the images which included the EVLA antennas was better by 
the expected amount under the assumption that they were identical to the standard VLA antennas. 
We note that the main benefit of the retrofitted antennas at 1.4 GHz is their performance over 
wider bandwidths, which will be available for use when the EVLA's new correlator comes online. 
For the current project, we proceed treating the retrofitted EVLA antennas no differently than the 
regular VLA ones. 

The final image was created by assembling the mosaic of these 22 individual pointing images 
using the task FLATN. For this mosaic, two 4096 x 4096 images with 1" pixels were created, centered 
on Coma 1 and Coma 3 (12^59™45!2 +27°57'53'.'l and 12^57™28!5 +2 7° 08' 08'.' 5, respectively). Only 
data out to the 30% power point of the VLA 1.4 GHz primary beam were used, which corresponds 
to the central 19 facets of each flys-eye. Data for each pointing were variance weighted according 
to the VLA primary beam power pattern. The rms sensitivity of the final image is fairly uniform, 
as shown in Figure [T] which depicts contours of constant rms as evaluated over 5' square boxes 
centered on each location. The majority of the formal Coma 1 and Coma 3 areas have rms below 
29 /iJy (Figure Ej), which correponds to a 5a detection of a Conia clus ter galaxy forming stars at 0.1 



Mq yr ^ (assuming the conversion of lYun. Reddy. Sz Condonll200ll ). The deepest regions within 



Coma 1 have an rms sensitivity of 22 /iJy, while those in Coma 3 are slightly less sensitive (Figure 
[3]). This is largely on account of the sidelobe response of the bright extended radio galaxies in the 
center of Coma 1 (NGC4874 and NGC4869). 



2.3. Source Catalog 

With varying sensitivity across the full mosaic image, we opted to create signal-to-noise (S/N) 
maps for source detection. We first calculated the local rms for each pixel in the final mosaic 
images, evaluated over a 5' box using the task RMSD (refer to Figure [H smaller-sized boxes are 
overly effected by real extended sources such as NGC4869). We then divided each mosaic image 
by its corresponding rms image to obtain the S/N maps. Source detection was then accomplished 
via the task SAD, which was directed to identify sources with peaks greater than 4.5o" and fit them 
with Gaussians. SAD was run on six separate portions of the full mosaic images to insure that 
sources were not missed due to the task's limitations for the maximum number of sources it will 
evaluate in a given run. We also had SAD create residual images, the purpose of which will be 



described shortly. The procedure is largely the same as that described in lBiggs &: IvisonI (j2006l ) 



The SAD catalogs then served as input for creating the final catalog. Each source in the SAD 
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catalogs was examined in the mosaic images (not the S/N image) using the task JMFIT, and 
its position, peak flux density, integral flux density and error, whether the source was resolved, 
and whether it was near the edge of the mosaic (within ~2f5, where the local rms measurement 
would be based on fewer resolution elements) were recorded. Repeated astrometric measurements 
with the VLA at 1.4 GHz indicate a positional accuracy of around C'l for strong unresolved point 
sources. The fainter sources in our survey have positional errors on the order of (/.'S (Equation 1 of 



Condon et al.lll998l ). Because JMFIT uses the noise measurement made over the full image when 
reporting flux density errors, we used the measured value at the source position in the RMSD map 
as the error in the peak flux density measurement. Similarly, we updated the integral flux density 
errors by scaling the JMFIT-reported values by the ratio of the local noise to the full mosaic area 
noise. Assessing whether a source is resolved is not entirely straightforward due to the mosaic 
nature of the images. Although greatly reduced by the smaller channels, bandwidth smearing does 
smear sources radially - and hence in different orientations for each pointing's contribution to a 
given source. Thus, an intrinsically unresolved source might appear to be marginally resolved on 
account of bandwidth smearing arising in the individual pointings that contributed to that source 
within the mosaic. We call a source resolved if the minimum size of its reported major axis from 
JMFIT was non-zero, however we note that there will exist false positives in the catalog. These are 
usually brighter sources, and can often be spotted through comparison of their peak and integral 
flux densities, which will be equivalent within their respective errors for unresolved sources. 

The residual images were then inspected to identify sources either missed or poorly fit by 
SAD. This was greatly simplified by the procedure of having run source detection on the S/N map, 
since contours at specific values of a could be plotted and thereby render the identification of such 
sources easy. Those sources missed by SAD were added to the final catalog after repeating the 
procedure outlined in the preceding paragraph. Because both SAD and JMFIT rely on Gaussian 
fitting, plotted a contours on the residual map also ease identification of sources for which fitting 
an individual Gaussian was problematic. In such cases, JMFIT was used to deblend sources which 
were better described as multiple overlapping Gaussians. Finally, for the more irregularly shaped 
extended sources (e.g., NGC4869) the task TVSTAT was used to perform aperture photometry. 
The catalog positions for these sources correspond to the location of their peak flux density. 

The resulting radio source catalog is presented in Table [2j There are over 1200 sources listed, 
although this figure includes lower significance detections and multiple entries for some sources 
which appear to be associated with a single extragalactic object (e.g., one entry for each lobe of 
a radio double). Not all sources that were fit by multiple components by SAD consist of multiple 
entries in Table [21 however. In cases such as NGC4874 where we have prior knowledge of the 
source and its origin in a single radio galaxy we have included only a single entry in Table [2j We 
have attempted to provide guidance for the possible multi-component sources by noting where our 
visual inspection of the images suggested a relationship between separate entries in the radio source 
catalog. This may be found in the "Comments" field of Table [2j The table includes sources down 
to S/N > 4.5(7, where S/N is based on the peak flux density and error from the SAD run on the 
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S/N map. Although 5a is traditionally used in radio surveys, we include sources with slightly 
lower significance. There are two reasons for this: First, they still have a very high chance of being 
non-spurious, especially in cases where they coincide with detected sources at other wavelengths. 
Second, diffuse sources with peak fiux densities less than five times the local noise might have 
integral fiux densities greater than five times their integral flux density errors. Restricting our 
acceptance criteria to sources with S/N > 5a and within the 0.92 square degree contiguous area 
formed by the Coma 1 and Coma 3 fields yields 628 radio sources. 



3. Matching to Optical Data 



The Sloan Digital S ky Survey (SPSS; York et^al. 20001 ) included the photometry for Coma in 
its Data Release 5 (DR5: lAdelman-McCarthv et a l. 2007). The SDSS observes the sky in drift-scan 
mode, using a novel multi-CCD camera (jCunn et al.jil998i ) that simultaneously collects data in five 
filters {u, g, r, i, and z). The median see ing at r is I'.'A, and the 95% com pleteness limit for point 
source detectio n is 22.2 mag in thi s filter (jAdelman-McCarthv et al.ll2007l ). The magni tudes are on 
the A B system ( Oke Gunn|[l983 ) and have uncertainties in the 0.02 - 0.03 mag range ( Ivezic et al. 
2004 ). Catalog data are a vailable online (http://www.sdss.org/dr5/[ ) and are described in the 



Early Data Release paper (jStoughton et al.ll2002l ). For correlation with the radio data, we used the 
SDSS DR5 catalog of all objects within large boxes that extended > 10' past each side of the formal 
Coma 1 and Coma 3 regions. These included both photometrically-classified stars and galaxies. 

We adopted a matching radius of 2" for sources with r < 22 and 1" for fainter sources. These 
figures were determined based on direct estimation using the data. First, Table [2] was correlated 
with the SDSS optical catalog, and the SDSS object nearest to each radio source was recorded. To 
evaluate whether a given radio-optical pair represented a real physical association, we repeated this 
correlation a number of times after having applied arbitrary fixed offsets to the radio catalog, thus 
arriving at estimates of the probability of chance superpositions derived directly from the data. The 
resulting distribution of radio-optical separations is shown in the top-left panel of Figure [H The 
large peak in the distribution of the unshifted data with small radio-optical separations is obvious, 
and the vast majority of such objects are real associations of radio emission with an optical source. 
At larger separations the two distributions are approximately the same, indicating that sources at 
such separations are consistent with expectations for random superposition of a radio source with an 
optical source. The false match rate appears higher than the actual data for most of the separations 
plotted in the top left panel of Figure [4] because the number of sources used in generating the shifted 
catalogs is conserved and the fraction of radio sources with real optical counterparts is high. In 
other words, if we were to remove the real counterparts and renormalize the unshifted data to have 
the same number of sources as the shifted data, the two curves would be more coincident. The 
two curves cross at a separation near 2", indicating that radio sources with optical counterparts 
closer than this value have a high probability of being physically associated with one another while 
radio sources with more distant optical counterparts are likely false associations. The density of 
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sources in the radio and optical catalogs is sufficiently high that at the larger plotted separations 
the number of matches declines. 

The use of the 1" search radius for fainter objects was determined in the same manner. Figure 
Hlalso includes the real and offset catalog results for three separate magnitude ranges, r < 20 (top 
right), 20 < r < 22 (bottom left), and r > 22 (bottom right). It can be seen that for each of 
the r < 20 and 20 < r < 22 samples the distributions of the real data and the shifted data cross 
somewhere just below 3". However, for the fainter sources (r > 22) the crossover is at separations 
less than 2". Figure S] also provides a direct way to estimate the confidence with which a given 
radio-optical association may be judged. For example, a pair with a 1" separation would have 
about a 1 in 40 chance of being a false association if the optical source had r < 20, about a 1 in 10 
chance if 20 < r < 22, and about a 1 in 5 chance if r > 22. A comparison of the three panels also 
shows the expected result that the optically brighter galaxies are more likely to be associated with 
a radio source. 

The S/N map that was used in the generation of the radio source catalog was then overlaid 
on the r-band images from the SDSS. The entire area of the radio mosaic was thus inspected to 
identify apparent matches of radio emission with optical sources, guided by an input listing of all 
SDSS sources with radio counterparts from Table [2] within 3". As indicated, we generally accepted 
sources with radio-optical separations up to 2" for brighter objects (r < 22, or Mr = —13 in Coma) 
but were more restrictive for fainter objects and used 1" separation as the cut. Some sources with 
separations smaller than these values were rejected for a variety of reasons. In some cases, multiple 
SDSS counterparts of a radio source existed and the radio emission was ascribed to the object 
that was deemed the more likely counterpart, with such a determination usually based on giving 
preference to the nearer counterpart with an SDSS object type of galaxy. Seventeen sources with 
separations greater than 2" were also accepted in cases where the radio emission was extended and 
its morphology was consistent with expectations based on the optical source. Examples include 
spiral galaxies with extended radio emission matching their optical disks (e.g., NGC4911) and 
eUiptical galaxies hosting "FRl" and "FR2" type radio sources (e.g., NGC4869, NGC4874). All of 
these exceptions are depicted in Figure O Finally, the visual inspection found four radio doubles 
where the midpoint of the two radio positions was within l" of the apparent optical counterpart. 

The full radio-optical source catalog is presented in Tabled and consists of 499 optical counter- 
parts to radio sources from Table [2j For completeness, we present the rejected optical sources within 
3" of a radio source in Table S] along with brief explanations for their exclusion from Table [3l The 
table includes summary data for the SDSS optical counterpart, including its photometric classifica- 
tion (galaxy or star) and ugriz "model" magnitudes. These are based on a common aperture defined 
by the more likely of a deVaucouleurs or pure exponential profile fit t o the data at r band, and 



are t he recommended magnitudes for color studies of galaxies (e.g., see lAdelman-McCarthy et al 



20071 ^ 
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4. Discussion and Summary 

Unless otherwise noted, in tiie ensuing discussion we restrict our analysis to the contiguous 
~ 0.9 square degree area defined by the combined Coma 1 and Coma 3 fields (see Section [1] and 
Figured]) and refer to this generically as the "survey area." 

About half of the radio sources within the survey area have optical counterparts in the SDSS. 
Of the 628 radio sources with S/N > 5a, 312 have optical counterparts from the SDSS catalog and 
an additional three optical counterparts to radio sources were found in our visual inspection of the 
radio data overlaid onto the SDSS r-band images (refer to Table O the optical sources not found 
within the SDSS catalog were presumably excluded due to their location near very bright objects 
such as diffraction spikes around saturated stars). A further 22 radio sources detected at less than 
5(T were also identified with SDSS optical counterparts within the survey area. 

Conversely, the fraction of optical sources (Ntot) detected by the radio observations (Njig) is 
shown in Figure[6l Specifically, we use the SDSS objects with photometric classification as galaxie^ 
and investigate the fraction of sources with > 5cr radio detections as a function of r magnitude. At 
the brighter magnitudes (r < 14, or < —21 for Coma members), our survey detects over half 
of the galaxies (15/28, including all six objects with r < 13.04). Velocity measurements for the 
radio detections, including many new measurements for fainter galaxies deter mined from an MMT 



campaign using Hectospec (Marzke et al. in preparation), are presented in iMiller et al.l (120081 ) 
where we derive and discuss the radio LF. We note here that NED contains measurements for most 
of the galaxies in the present discussion (those with r < 17.5) since Coma has been the targe t 



of numerous galaxy spectroscopic surveys, including the SDSS (lAdelman-McCarthv et al.l l2008l) 



Assu ming a redshift range 4,000 km s ^< cz < 10,000 km s ^ for Coma (e.g., IColless Sz Dunn 



19961 ). all 15 of the radio detections with r < 14 are found to be Coma members. The radio 
detection fraction remains in about the 10% - 20% range for 14 < r < 19 before dropping. At 
brighter optical magnitudes the radio sources are almost exclusively Coma members, while at the 
fainter magnitudes they are predominantly background galaxies. By r ~ 21 the radio-detected 
fraction falls to about 1%. 

The types of the radio-detected sources are elucidated by their optical colors. Figure [7] shows 
the optical color-magnitude diagram, with radio detections indicated by filled circles for cluster 
members, open circles for non-members, and open triangles for sources without velocity information. 
The red seq uence is very clear l y del ineated, and we fitted it for galaxies with r < 17 using the 



procedure of lLopez-Cruz et al.l (j2004l ). Note that this fitting is strictly based on photometric data, 
so some galaxies that lie on the red sequence will actually be non-cluster members whose colors 
happen to be consistent with the Coma red sequence. It can be seen from Figure [7] that the radio 
detection fraction for red sequence galaxies rises at r = 14.5 {Mr = —20.5 for Coma). Setting the 



^Six bright objects (r < 12) were removed from the optical catalog as inspection of the images indicated they were 
saturated stars mistakenly classified as galaxies. 
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red sequence range at ±2(7 about the fitted relation of {u — g) = 3.192 — 0.095r, 18 of 43 red sequence 
galaxies with r < 14.5 are radio detect ions. The detection ra te increases with optical luminosity, as 
found in numerous other studies (e.g.. iLedlow &: Owenlll996l ). All five red sequence galaxies having 
r < 13.2 {Mr < —21.8) are detected by the radio data (these include four of the five brightest cluster 
galaxies noted above). Conversely, only two of the 135 red sequence galaxies with 14.5 < r < 17 
(—18 > Mr > —20.5) are radio detections. Based on velocities from the literature, one of the 
18 radio detections among the brighter galaxies is a background source (J12581169+2708117 with 



cz = 11889 km s : lAdelman-McCarthy et al.ll2008l ) alon g with one of the two ra dio-detected fainter 
galaxies (J12581684+2710267 with cz = 18179 km s~^: IColless DunnI figgg ) . This means that 
the faintest radio-detected Coma niember on the red se quence has Mr = —19.9, with the possible 
exception of J12591536+2746052. iBiviano et al.l ( 19951 ) report a velocity of 6315 km s~^ for this 
galaxy, although they indicate this determination is of intermediate quality as it is based on few 
spectral features. If true, it would have Mr = —14.8 and L1.4 = 2.1 x 10^'^ W Hz~^, making it an 
extreme outlier about four magnitudes fainter than the next faintest Coma radio-detected member. 

The detection fraction is similar for galaxies bluer than the red sequence. Of the 66 galaxies 
with r < 17 and (u — g) colors at least 2a bluer than the fitted red sequence, 25 are radio detections. 
Most of the detections have r < 16, where 17 radio detections are found among 34 total galaxies. 
Again, relying on published velocity measurements from NED and the SDSS we find that 6 of 
the 25 radio-detected galaxies are background galaxies (2 of 17 for the r < 16 detections). The 
faintest radio-detected Coma member with a blue color has Mr = —18.8. Presumably, the radio 
emission from these blue Coma galaxies has its origins in star formatio n. This is lar g ely co nfirmed 
via further analysis of photometric data and optical spectroscopy in iMiller et al.l ([200a), where 
only one blue Coma galaxy showed spectroscopic evidence for an AGN in the form of emission 
line diagnostics typical of a Seyfert galaxy (ARK395). Its blue colors (including UV) extending 
across its disk suggest that star formation also contributes to the net radio emission from this 
galaxy. Based on their respective radio and optical morphologies, some of the blue Coma galaxies 
appear to be undergoing ram pressure stripping. For example. Figure [5] includes IC4040 which has 
radio contours compressed to its northwes t side and extended to th e southeast - consistent with 
the H I deficit and morphology reported in iBravo-Alfaro et al.l (|200d ). This disturbed morphology 
produces an offset between the fitted radio position and the nucleus of the galaxy and hence its 
inclusion in Figure O Although it lies outside the Coma 3 footprint, KUG1255+283 (also shown 
in Figure [5]) has a similar radio morphology, in this case with its radio emission extended away 
from the galaxy in a direction roughly opposite that of the cluster center. Other examples of known 
star- forming galaxies with radio contours suggestive of ram pressure interactions include Markarian 
57 and KUG1258+279A with radio-optical positional offsets of 1'.'3 and I'.'O, respectively. Again, 
in each case the radio contours are asymmetric with respect to the galaxy nucleus and extended 
in a direction roughly po inting away from the center of the cluster. Both are included in the 
Bravo-Alfaro et al.l pOOol ) work, although Markarian 57 is actually H i rich. KUG1258+279A is 
H I deficient, although o nly weakly detected and hence its H i morphology is only suggestive of 



ram pressure interaction (IBravo-Alfaro et al. 



2001 
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Although this survey specifically addresses the Coma cluster, its combination of sensitivity 
and areal coverage are comparable to some published deep fields. The total area covered within the 
combined Coma 1 and Coma 3 footprint is about 0.9 square degrees, with a typical rms sensitivity 
of 28 /iJy per A'.' A beam (refer to Figures □ and ED . The VLA-VIRMOS Deep Field (VVDS) includes 
radio observa.tions with an rms sensitivity of 17 /iJy per 6" beam across a 1 square degree region 



(jBondi et al.l 120031 ) . and the Subaru/XMM-Newton Deep Field (SXDF) covers 0.8 square degrees 



with 20 /iJy sensitivity at about 5" by 4" resolution (note that these figures are specifically for 



the s ubset of the total area surveyed having local rms sensitivity below 20 //Jy; ISimpson et al 



20061 ). Along with the excellent multiwavelength data available for Coma (refer to Section [T]) 
including optical spectroscopy of galaxies down to fairly faint magnitudes (r ~ 21; Marzke et al., 
in preparation), the multi- waveband Coma data will also be of value to deep field science. 

Figure M shows the source counts for our survey, and the source count data are provided in 
Table m We have restricted this analysis to the 628 sources within the Coma 1 and Coma 3 coverage 
area that were detected with a peak flux density greater than five times the local noise. In addition 
to source counts fro m the VVDS and SXDF, we also plot those of the Phoenix Deep Field (PDF; 



Hopkins et al.ll2003l ) and the FIRST survey. The wider area of these surveys improves the numbers 
and hence error bars for source counts at higher fiux densities. We note that the source counts 
for the Coma data have only been corrected for the areal coverage in each bin (i.e., sources in 
the fainter bins are detectable only in the deeper regions of the survey area), thus ignoring effects 
such as Eddington bias and resolved sources with integral flux densities sufficient for inclusion but 
peak fiux densities failing to reach the 5a cutoff. Even so, the source counts are consistent with 
these other studies, showing the usual decline in normalized counts between 100 mJy and 1 mJy, 
followed by leveling below 1 mJy. A slight excess of sources is seen in the 1 to ~ 5 mJy region, 
and it is tempting to attribute this to Coma member galaxies since these flux densities correspond 
to luminosities of a few times 10^^ W Hz^^ for galaxies at the distance of Coma. While Coma 
galaxies certainly do contribute to the source counts between ~ 1 to ~ 5 mJy, the excess over the 
FIRST survey source counts in this region is much greater than can be explained by the presence 
of the Coma cluster. For the source counts at approximately 1.5 mJy and 3 mJy in Figure [HI 
the difference between the number of sources in the Coma survey area and that expected based 
on the FIRST source counts amounts to roughly 20 to 25 sources per each bin. The matching to 
SDSS optical data and optical spectroscopic databases indicates that only nine Coma members have 
radio flux densities falling in these bins (three for the 1.5 mJy bin and six for t he 3 mJy bin ) . The 



comprehensive optical spectroscopic coverage for radio-detected galaxies (see iMiller et al.l l2008l ) 



insures that we are not missing Coma galaxies down to very faint optical magnitudes {Mr < —15). 
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Fig. 1. — Contours of constant rms sensitivity for final mosaic image. Rectangular areas for Coma 
1 (upper left) and Coma 3 (lower right) are indicated, along with crosses marking the locations of 
the 22 VLA pointing centers. The innermost contour for each Coma 1 and Coma 3 is at 25 ^3y, 
with successive contours at 30, 35, 40, 45, 50, and 55 ^Jy. 



Fig. 2. — The (u,v) coverage for the central pointing in the Coma 1 mosaic (Comal_F in Table [T]). 
Every 51** record is plotted. 
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Fig. 3. — Fractional area covered at a given sensitivity limit for Coma 1 (dotted line), Coma 
3 (dashed line), and combined (solid line). The top axis is the associated star formation rate 
for a source detected at Sa, d etermined using the 1.4 GHz star formation rate calibration from 



Yun. Reddv. k CondonI ((20011). 
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Fig. 4. — Number of radio-optical pairs as a function of separation. The results for the actual radio 
catalog are plotted as filled circles, while the frequency of chance superpositions as evaluated by 
applying shifts to the radio catalog are shown as open circles. Top left: for all SDSS sources; top 
right: for SDSS sources with r < 20; bottom left: for SDSS sources with 20 < r < 22; bottom 
right: for SDSS sources with r > 22. 



Fig. 5. — Radio galaxies with separations over 2". Each square image is 2' by 2', with radio 
contours at 3, 5, 8, 13, ... times the local noise. The top right galaxy is IC3913, and the middle 
left is NGC4839. The galaxy at the bottom left is KUG1255-I-283, and appears to be undergoing 
ram pressure stripping. 
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Fig. 5. — continued. The galaxy in the middle right panel is NGC4874, and the bottom panel is 
NGC4869. 

Fig. 5. — continued. The galaxy in the middle right panel is IC4040, which appears to be undergoing 
ram pressure stripping. The bottom left galaxy is NGC4911. 
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Fig. 6. — Histogram of r magnitudes for SDSS galaxies with radio counterparts for sources within 

the Coma 1 and Coma 3 area (left vertical axis indicates number of radio galaxies, A^rg)) and 
fraction of total SDSS galaxies {Nrq / Ntot) detected by radio observations within that magnitude 
bin (right vertical axis). 
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Fig. 7. — Color magnitude diagram for all SDSS galaxies within the Coma 1 and Coma 3 regions. 
Non-radio detected galaxies are shown as small grey crosses, while radio-detected galaxies are 
indicated by larger symbols: filled circles are spectroscopically-confirmed Coma members (4,000 
km s-i< cz < 10,000 km s"^), open circles are spectroscopically-confirmed background galaxies 
(with the exception of the blue object at r = 15.5, which is the planet ary nebula PG1257 +279), and 
open triangles are sources for which no spectroscopy is available (see lMiller et al.l (|2008l ) for details 
on spectroscopy of Coma radio-detected galaxies). The long dashed lines are the ±20" boundaries 
to the fitted Coma red sequence, while the vertical dotted line shows the approximate limit for 
radio detection of red sequence members. 
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Fig. 8. — Source counts for Coma data (filled black triangles) along with those derived from similar 
surveys. 
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Table 1. Observation Summary 



Pointing ID RA(J2000) Dec(J2000) Dates Observed 
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Note. — The pointings for each field are arranged in three north-south columns 
where A-D-G-J form the first column, C-F-I from the middle column, and B-E-H-K 
form the last column. Th e "F" pointings thus cor respond to the formal Coma 1 and 
Coma 3 field centers from iKomivama et al.l (j2002l ). 
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Table 2. Radio Source Catalog 



ID 


RA(J2000) 


Dec(J2000) 


S/N 


Speak 


rms 


Sint 




Res? 


Edge? Comment 


C3A-001 


12 55 11.22 


27 27 22.4 


7.9 


921 


115 


1880 


313 


1 


1 


C3A-002 


12 55 12.38 


27 27 39.8 


4.8 


556 


114 


570 


201 


1 


1 


C3B-917 


12 55 12.46 


26 59 39.6 


4.6 


551 


119 


1082 


316 


1 


1 


C3B-001 


12 55 12.57 


26 58 58.4 


12.6 


1516 


120 


3562 


443 


1" 


1 


C3B-003 


12 55 13.30 


26 59 48.4 


8.5 


1007 


118 


2236 


358 




1 



''Extended source with integral flux density and error determined in irregularly-shaped aperture using 
TVSTAT. 

Note. — ONLY A PORTION OF THE TABLE IS SHOWN HERE. Column descriptions: (1) unique 
source ID that we found useful in our analyis; (2) & (3) right ascension and declination (J2000) of radio 
source fitted peak; (4) S/N of fitted peak as determined by source identification algorithm (SAD) on S/N 
map (see text). Sources with 4.94 < S/N < 5.0 are listed as S/N = 4.9 to clarify that they are not formally 
5(7 detections.; (5) peak flux density in fiJy; (6) local rms as evaluated in a 5' square box centered on the 
source position (this is also the error associated with the peak flux density measurement); (7) & (8) integral 
flux density and its associated error, each in /iJy; (9) binary code indicating whether the source was resolved 
{— 1) or unresolved {— 0); (10) binary code indicting whether the source was within ~2.'5 of the edge of the 
radio mosaic (= 1); and (11) comments. 
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Table 3. Optical Counterparts to Radio Sources 



RA 


Dec 


Sep 


u 


9 


r 


i 


z 


Class 


Radio ID 


Comments 


12 55 17.18 


26 45 48.2 


0.39 


24.956 


24.918 


23.207 


22.029 


21.345 


6 


C3C-002 




12 55 20.64 


27 14 09.2 


0.72 


22.006 


21.768 


21.761 


21.148 


21.238 


6 


C3B-006 


C3B-005, C3B-009? 


12 55 23.33 


27 02 13.9 


0.48 


23.073 


21.895 


20.626 


20.015 


19.275 


3 


C3B-011 




12 55 23.90 


26 56 40.9 


0.71 


24.052 


21.926 


20.222 


19.377 


18.987 


3 


C3C-003 




12 55 24.10 


26 57 37.4 


0.41 


23.894 


24.526 


24.003 


22.211 


21.630 


6 


C3B-012 





Note. — ONLY A PORTION OF THE TABLE IS SHOWN HERE. Column descriptions: (1) Right ascension of SDSS 
optical counterpart (J2000); (2) Declination of SDSS optical counterpart (J2000); (3) Separation between radio position and 
SDSS optical counterpart, in arcsec; (4 - 8) SDSS ugriz model magnitudes; (9) SDSS photometric class of optical counterpart, 
where 3 = galaxy and 6 = star; (10) Radio ID, from Table [51 (11) Comments, where: 

a = Radio emission is resolved, providing a likely explanation for the somewhat large radio-optical separation, 
b = The morphology of the extended radio emission is consistent with that of the optical. In the case of spiral galaxies this 
indicates that the radio emission traces the galaxy disk, whereas for elliptical galaxies the center of the galaxy appears to be 
where the radio emission originates. 

c = Source appears to be the host of a radio double. The indicated separation is for distance from optical source to midpoint 
of the radio counterparts listed in columns 10 and 11. 
d = Source separation is greater than 3". See Figure [5] 
e = Optical position for bright clump within larger galaxy, 
f = Galaxy appears to be one member of a galaxy pair. 

g = Radio emission is unresolved, but lies within the optical extent of the galaxy. 

h = Optical source not in SDSS catalog due to proximity to bright source, usually the diffraction spike of a saturated star. 
The position represents the coordinates of the peak of the object as measured directly in the SDSS r band image. 
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Table 4. Rejected Optical Counterparts to Radio Sources 



RA(J2000) 


Dec(J2000) 


Sep 


u 


9 


r 


i 


2 


Class 


Radio ID 


Comments 


12 55 33.15 


26 41 14.3 


1.04 


23.394 


22.900 


21.996 


21.205 


20.494 


3 


C3C-007 


a 


12 55 45.55 


26 50 05.3 


1.27 


26.666 


23.608 


21.443 


20.295 


19.620 


3 


C3C-010 


a 


12 55 47.06 


27 27 50.8 


2.32 


21.903 


20.786 


19.543 


19.016 


18.544 


3 


C3A-013 




12 55 49.34 


27 33 30.6 


2.49 


23.872 


21.690 


20.717 


20.523 


20.023 


3 


C3A-015 


b 


12 55 58.32 


27 21 23.0 


2.45 


24.888 


24.555 


22.777 


21.511 


20.637 


6 


C3A-020 


b 



Note. — ONLY A PORTION OF THE TABLE IS SHOWN HERE. Column descriptions: (1) Right ascension 
of rejected SDSS optical counterpart (J2000); (2) Declination of rejected SDSS optical counterpart (J2000); (3) 
Separation between radio position and SDSS optical counterpart, in arcsec; (4 - 8) SDSS ugriz "cmodel" magnitudes; 
(9) SDSS photometric class of optical counterpart, where 3 = galaxy and 6 = star; (10) Radio ID, from Table[21 (11) 
Comments, where: 

a = Although the radio-optical separation is acceptable for the r magnitude of the potential counterpart, visual 
inspection of the overlaid radio contours strongly suggested the sources are unassociated. 
b = The radio source was assigned to a galaxy with a smaller separation (see Table 
c = Same comment as "a" plus the radio source is detected at less than 5o-. 

d = The radio morphology is that of a radio double source, and this optical counterpart appears to be a chance 
superposition of a galaxy with one of the lobes. The optical counterpart for the radio emission would be expected to 
lie in between the lobes, and not coincident with one. 

e = Visually appears offset. The radio emission also appears to be a lobe of a separate FR2 type source, and hence 
not associated with this galaxy. 

f = This radio source was assigned to a brighter galaxy with radio-optical separation greater than 3". It may be seen 
in Figure [S] 



Table 5. Source Counts within Coma 1 and Coma 3 







s 


< 5 > 


N 


Factor 


S'^-^dN/dS 




[mJy] 


[mJy] 






[Jyi-^sr-i] 


0, 


,110 


- 0.140 


0.129 


50 


2.81 


3.20± 0.45 


0, 


,140 


- 0.180 


0.158 


109 


1.32 


4.04± 0.39 


0, 


,180 


- 0.230 


0.202 


87 


1.02 


3.71± 0.40 


0, 


,230 


- 0.280 


0.253 


73 


1.00 


5.32± 0.62 


0, 


,280 


- 0.370 


0.323 


71 


1.00 


5.30± 0.63 


0, 


,370 


- 0.650 


0.483 


75 


1.00 


4.91± 0.57 


0, 


,650 


- 1.200 


0.889 


38 


1.00 


5.83± 0.95 


1, 


,200 


- 2.000 


1.546 


40 


1.00 


16.8 ± 2.7 


2, 


,000 


- 4.000 


2.969 


43 


1.00 


37.0 ± 5.6 


4, 


,000 


- 10.000 


5.757 


24 


1.00 


36.0 ± 7.3 


10, 


,000 


- 99.999 


35.783 


15 


1.00 


145. ±37. 
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